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FOREWORD 


This is the Final Report on Contract NAS9-12055 which 
was executed by Vought Systems Division between January 1972 
and October 197^. This is a Summary Report; the information 
contained in this report has been previously presented in 
detail in Progress Reports, and in two technical papers which 
are included as Appendices to this report . 

The NASA-JSC technical monitors on this program were 
Mr. D. W. Morris on the system demonstration phase, and Mr. 
Burrell French on the zero-g condensation experiment phase. 
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In 1970 and 1971 Vought Systems Division (VSD) performed a study on 
use of refrigeration in rejection of future spacecraft environmental control 
system heat loads, under NASA-JSC contract NAS9-9912. One conclusion of that 
study was that vapor compression was the most attractive mechanical refrigeration 
concept for use in future spacecraft (Reference [l] ) - A revolutionary heat 
rejection system of a hybrid type was recommended for further study (Reference [2]). 
This was a hybrid refrigeration/radiator system which acts as either a refrigeration 
system or as a radiator system, depending on the vehicle heat rejection demand 
and the external environment. The hybrid system uses a common flow network 
and radiator condenser; compressors and pumps are valved into and from the flow 
loop as required. This system is particularly well suited to lunar surface shelters. 
Peak heat rejection requirements during the long lunar day coincide with peak 
power availability from solar cells; during the lunar night while power is scarce, 
the relatively efficient radiator mode is used. 

There were several technical questions relating to the feasibility of this 
hybrid concept, including: 

(1) Could the system be changed from the single-phase "radiator" 
mode of operation to the two-phase "refrigeration" mode? 

(2) Could a radiator /condenser design be evolved which would 
provide satisfactory operation in both the "radiator" 
and the "refrigeration" modes of operation while in a 
zero-g field? 

(3) Could a highly efficient motor-compressor unit capable 
of operation in a zero-gravity field be obtained? 

The first of these questions was auickly answered in a feasibility demon- 
stration test (Reference [3]), where it was shown that a hybrid system could be 
switched from one operating mode to the other without damage to system components. 

NASA-JSC then contracted with VSD (in the subject contract) to design 
a suitable radiator/condenser and to demonstrate it in a thermal-vacuum test in 
a systems test. Later this contract was modified to include design, fabrication 
and ground checkout of an experiment apparatus which would verify the condenser 
flow design under simulated zero-g conditions. 

This report briefly describes the design, fabrication and testing which 
were accomplished under contract NAS9-12055. 
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DISCUSSION 


Systems Test 


The purpose of the systems test was to demonstrate that a radiator/ 
condenser could be designed which was efficient in both condensation ("refrigera- 
tion") and liquid phase ("radiator") operating modes. In addition, switchover 
from the "refrigeration" mode to the "radiator" mode and vice versa was to be 
demonstrated. 

The major issues associated with design of the radiator/condenser 


(l) tube layout 
and (2) tube diameter 

There was a desire to use the wide heat load range radiator concept (References 
[U], (5l» and [6]) which has tubes of varying length as shown in Figure 1. In 
a liquid phase radiator this approach provides a heat load ratio of 100 to 1 
or more, which is a distinct advantage in a spacecraft heat rejection system. 
However, it was found that the unequal flow lengths of the wide heat load 
radiator were unsuitable for a condensing radiator. There were two primary 
reasons for this: 

(1) Orif icing required in the wide heat load radiator tubes 

is difficult to size so that proper flowrates are established 
in both liquid and vapor flow. 

(2) The wide heat load radiator has inherently higher pressure 
drop than a conventional parallel flow radiator, and this 
degrades performance in the two-phase or "refrigeration" 
mode of operation. 

It was found that suitable performance could be obtained in both 
"radiator" and "refrigeration" modes of operation with tubes of 3.175 mm (l/8-in) 
inside diameter. This diameter tube, with the appropriate flowrates in each 
operating mode, produces acceptable pressure drop and heat transfer rates. In 
addition, it also promised flow stability in the condensing mode of operation 
in zero-gravity. In order to accomplish the design of the condensing radiator 
it was necessary to predict the pressure drop of a condensing two-phase flow 
in zero-gravity. A method for doing this was developed, based on the Lockhart- 
Martinelli approach (References [7] and [8]). The method developed is discussed 
in detail in Reference [9] and in Appendix A. It vis also necessary to predict 
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FIGURE I WIDE HEAT LOAD RANGE RADIATOR TUBE LAYOUT 


3 



the flow regime which would prevail along the condensation path. The 
method developed by Baker (Reference [10]) was used to do this as explained 
in Appendix A. 

Two condensing radiator panels were fabricated as shown in Figure 2, 
and as described in Appendix A and in Reference [ll]. A dry expansion evaporator 
was designed and built as shown in Figure 3- The hybrid refrigeration system 
was assembled with the two radiator/condenser panels installed in the VSD 
Space Environment Simulator. A schematic of the system is shown in Figure U, 
and photographs of the installation are given in Figures 5-8* Testing 
of the system was accomplished under both atmospheric and vacuum conditions as 
described in Reference [8], Typical operating conditions are given in Tables 1 
and 2. 

Results of the testing were entirely satisfactory. Typical temperature 
profiles along the length of the tube are shown in Figures 9-12. figures 9 
and 10 clearly show the break-point in the temperature profile which occurs at 
the point of complete condensation. While a two-phase condition is maintained 
the fluid temperature is a function of the saturation pressure, and falls 
slightly in the direction of flow due to flow pressure drop. When a single- 
phase flow is established, then the temperature is independent of pressure, 
and so falls in proportion to the panel heat transfer rate. The predicted 
temperatures compare favorably with the data; the abrupt temperature changes 
in the predictions near the entrance are a result of longitudinal conduction 
not being considered in the mathematical model. The large temperature drop 
results from the small heat transfer coefficient near the end of the super- 
heated gas flow section of the condenser. As condensation starts in the two- 
phase flow region, heat transfer to the tube increases dramatically. 

Pressure drop predictions were generally higher than the data with 
an error of generally less than 20 % as shown in Table 3- This is good correlation 
for two-phase flow pressure drop (Reference [8]). 

There were no signs of flow instability in these tests. Runs were 
made with the radiator panels tilted and with vibration imposed cn the panel 
in attempts to induce flov instability. At no time was there any indication 
of such instabilities, however, the nature of the test hardware and instrumen- 
tation could obscure instabilities. 
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TABLE 1 

NOMINAL OPERATING CONDITIONS FOR REFRIGERATION SYSTEM MODE TESTS 


Refrigerator 

Refrigerant 

COP 

Condenser Temperature 
Evaporator Temperature 
Evaporator Heat Load 
Condenser Heat Rejection 
Compressor Gas Outlet Temperature 
Refrigerant Flowrate 
Condenser Sink Temperature 

Evaporator/ Intercooler Heat Load Flow Loop 
Heat Transport Fluid 
Inlet Temperature 
Outlet Temperature 
Flowrate 


R12 

1.75 

140°F (p = 221.32 psia) 
40°F (p = 51.67 psia) 
6500 BTU/HR 
11300 BTU/HR 
160°F 
180 LB/HR 
-300° F 


Water 

100°F 

87°F 

500 LB/HR 
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TABLE 2 

NOMINAL OPERATING CONDITIONS FOR PUMPED FLUID RADIATOR MODE TESTS 


e 

Radiator 


0 

Transport Fluid 

R12 

it 

Inlet Temperature 

'u 

o 

O 

u 

Outlet Temperature 

40° F 

II 

Flowrate 

400 LB/HR* 

L) 

Heat Rejection 

8000 BTU/HR 

u 

Sink Temperature 

-300° F 

i 

i ii_ 

Pump Outlet Pressure 

150 psia 

: Li 

Evaporator/Intercooler Heat Load Flow Loop 


: i i 
Li 

Transport Fluid 

Water 


Inlet Temperature 

180°F 

* i \ 

i 

Outlet Temperature 

164°F 

t 

i ! 

FI owrate 

500 LB/HR 


m 


u 

o- 

0 


* This yields R12 flow at 44.5 LB/HR/TUBE for 9 parallel tubes on a 0.063 Inch 
sheet with 8 Inch centers across a 6 ft span which Is essentially equivalent 
to the intended two-dimensional Modular Radiator System Test Condition for 
R21 flow through twelve tubes on a 0.020 Inch sheet with 6 Inch centers , 
Reference 5. 
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FIGURE 10 RC PANEL TEMPERATURE PROFILE 

(EVAPORATOR HEAT LOAD - 9250 BTU/HP.) 
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TABLE 3 

PRESSURE DROP CORRELATION 


PANEL 

NO. 

P IN 

t in 

W F-12 

4P cal 

ap test 

ERROR 

#1 

n 

PSIG 

°F 

lbm/hr 

PSI 

PSI 

% 

X 


167 

126 

21.4 

2.9 

2.8 

3.6 

X 


152 

118 

20.0 

2.9 

2.5 

16.0 

X 


172 

178 

18.2 

2.1 

1.7 

23.4 

X 


123 

152 

16.9 

2.1 

1.7 

23.4 

X 


205 

181 

23.1 

2.8 

2.4 

16.7 


JL 

237 

. .185. „ 

25.4 

11.5 

... .8.4 _ 

J7J1 


X 

240 

150 

26.7 

10.4 

8.9 

,1_6._8 


X 

168 

157 

20.0 

10.1 

7.0 

43.1 


. X_ ... 


152.. _ 

. 20.0 

..8,5 

8c_4 _] 




_X 

134 

132 

16.7 

_ .7 ._4 

. „7 t 2_ j 

...L.8 


PANEL #1 - TUBE I.D. » 0.125 in. 
PANEL #2 - TUBE I.D. * 0.095 in. 
AP CAL ‘ CALCULATED PRESSURE DROP 
AP test - ACTUAL PRESSURE DROP 
P JN - RC PANEL INLET PRESSURE 
T in - RC PANEL INLET TEMPERATURE 
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The overall performance of the system was within expected bounds. 

Figure 13 shows the Coefficient of Performance (COP) measured in the system 
as a function of compression ratio compared to prediction. There is excellent 
agreement here; at the low shaft speeds (200 rpm) the efficiency of the 
compressor falls considerably below optimum, as the data shows. The predic- 
tion is for optimum compressor speed at each load condition. It should be 
noted that this curve gives COP exclusive of the drive motor. Figure Ik 
shows a performance map for the system tested. 

The last objective of the test was to demonstrate that switchover 
from the single-phase ("radiator") mode of operation to the two-phase ("refrigera- 
tion") mode and vice versa could be made routinely without damage to the system. 
Numerous switchovers were made according to the operational procedures developed 
in the test plan, and no difficulty was encountered, nor was any damage to 
hardware sustained. 

Finally, a considerable amount of experience in operating the hybrid 
refrigeration system was obtained. 

The systems test is discussed in much greater detail in Reference [ 1 1 ] . 
2.2 Zero-G Experiment 

The system tests left unanswered the questions regarding heat transfer 
rate, pressure drop, and flow stability of condensing flow of a wetting fluid 
in a zero-gravity field. In order to increase confidence in the design, and 
to reduce technical risk, it was decided to design an experiment to verify the 
analyses in simulated zero-gravity. 

An investigation of the various means of achieving simulated zero- 
gravity was conducted, and it was decided that the most suitable approach was 
through the use of an aircraft flying parabolic trajectories. This approach 
can produce low acceleration conditions for periods of up to 30 seconds. 

The experimental approach which evolved was to film the condensing 
fluid flow with a high speed camera in local regions, and to examine the flow 
down the entire length of a tube visually. 

The design of the zero-g test apparatus is discussed in detail in 
Reference [12]. The system was designed for an operating pressure of 1.72 MPa 
(250 psia) at a temperature of 378°K (220°F), with delivery flowrates up to 
25 g/s (200 lb/hr). The system was tested operating !n «teady state at 
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FIGURE 1? COEFFICIENT OF PERFORMANCE VS PRESSURE RATIO 




pressures up to 1.&* MPa (238 psia), and was operated in a transient 
up to a ores sure of 2.07 MPa (300 psia) vhere pressure relief valve 
operation was verified. 

Three fused quartz tubes from the lot used in the flow visuali- 
zation apparatus were pressurized to failure; all burst at a pressure of 
greater than 13.8 MPa (2000 psia). The freon tank was fabricated to ASME 
unfired pressure vessel code, and was proofed to U.13 MPa (600 psi). The 
gas receiver was designed for a burst pressure of 6.89 MPa (1000 psi) and 
was proof-tested at 3.U5 MPa (500 psi). The structure of the test apparatus 
was designed for loads of l6 g's forward along the longitudinal axis, 5 6' 9 
down and 2 g’s in other directions. Figures 15-20 show the test apparatus 
in some detail. 

The tests conducted on the ze~o-g condensation experiment apparatus 
are described in detail in Reference 1X3] • and are also described in 
Appendix A. The flight tests are described in detail in Reference [lU] and 
in Appendix A Results of these tests are discussed in Appendices A and B. 
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FIGURE 17 CLOSE-UP VIEW OF PLEXIGLASS FLOW BOX 




















3.0 CONCLUSIONS AND RECOMMENDATIONS 

3.1 Conclusions 

The primary conclusions of this work are: 

(1) The hybrid conventional radiator /vapor cycle heat rejection 
system is feasible. 

(2) Operation of a spacecraft refrigeration system is feasible. 

(3) Flow regime (e.g., annular, slug, bubble, etc.) in condensing 
flow in zero-g can be predicted reasonably well with the Baker plot. 

(It) Condensing flow heat transfer rates and friction pressure drops 
do not appear to be degraded under zero-g conditions, and are reasonably well 
predicted using conventional techniques. 

(5) Condensing radiator can be designed so that hydrodynamic ally 
stable flow conditions prevail. 

3.2 Recommendations 

The following recommendations are made relative to this work: 

(1) Additional simulated zero-g experiments should be conducted 
under more carefully controlled conditions and with better instrumentation. 

(2) Data taken in simulated zero-g conditions should be used to 
improve techniques for calculating heat transfer rates and pressure drop. The 
photographic data taken in the subject tests clearly shows that the flow 
models used in conventional correlations are not correct. 

(3) Attention should be directed to the practical design of other 
system components for use in zero-gravity, including the dry-expansion evaporator 
and compressor lubrication. 
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APPENDIX A 


DEVELOPMENT OF A DIRECT CONDENSING RADIATOR FOR USE 
IN A SPACECRAFT VAPOR COMPRESSION REFRIGERATION SYSTIM 
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Development of a Direct Condensing 
Radiator for Use in a Spacecraft Vapor 
Compression Refrigeration System 

The development and feasdnlity itsting of a hyhrul spacecraft heat rejection system that 
incorporates a single raduitor captible of functioning as either u conventional space 
radiator or as a condenser in a refrigeration cycle is described Emphasis ts placed on 
development of the raduitor condenser Rid, which is considered to be the most critical 
component of the hybrid s • stem The selection, design and fabrication of candidate RC 
configurations are described together with preliminary parametric analyses necessary to 
establish pressure drop, heat transfer and How stability characteristics Verification 
testing in one-g ami zero-g environments ts desrcibed, the litter condition being obtained 
by means of a C115 aircraft The testing included How rts utilization 'i e , high speed 
photography ) of the condensation processes tu a parallel channel qiuirtz tube system 
modeling the RC Representative qualitative photographs are presented Results indi 
cate stable flow conditwns prevail for both one-g and zero-g operation 
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latniactiaa 

The objective of the praeent Mudy wa* to develop u 
condensing space radiator for use m a spacecraft vapor compres- 
sion refrigeration system. The present study wa» part of a pro- 
gram to develop a hybrid spacecraft heat rejection system with a 
single radiator that could function as either a conventional -pace 
radiator or as a condenser in a refrigeration cycle, thus the space 
radiator flow parages must be capable of accommodating either 
condensing two-phase flow or single-phase liquid flow 
The basic approach taken involved feasibility testing to demon- 
strate the operation of the hybrid heat rejection system, followed 
by development of individual component" The most critical 
component was considered to be the space radiator condenser 
(RC }, which is the subject of this paper 

Development of the KC included <1 > a requirement* study to 
define the likely configuration of advanced spacecraft radiators, 
(2) parametric analyses of RC’* to establish pressure drop, heat 
transfer, and flow stability characteristic-, >A\ -election and de- 
sign of candidate configuration", <4 fabrication of candidate 
configurations, i5) design and fabrication <>f an experiment 
package for cinematography of condensing two-phase flow, and 
tti) verification testing in one-gravity and m Jtero-g whib fly ing 
Keplerian trajectories in a f-135 aircraft 

lulgraiil 

Refrigeration of the Environmental Control >v-tem M> 

Contributed b\ the Aerospace Division and presented at the Inter- 
•octoty Conference on hnur»n mental sterna San Die**., < ahf . 
Ju)> 1ft 19, 1973. of The Aulnic \n Sck it n o» Mkhwou In 
oiNECBa Manuscript received at AsMIl Headquarters April 17 
1973. Paper No 73-K\ Aa-5 


heat load in manned spacecraft has not been neceasary to date 
Kxpendable heat sink* such as water evaporators and conven- 
tional pumped fluid radiator* have t>een entirely adequate 
However, the following factor-, either individually or collec- 
tively, may alter the presently satisfactory .ituation 

1 Construction of large space stations and other structure*, 
coupled with clu"termg of many vehicle- together t.such aa in 
docking a resupplv vehicle into a space -tation i may dramatically 
reduce the accessibility of the excellent -pace heat sink at lime* 

2 Internal power requirements, and thu.- heat rejection re- 
quirements, may grow at a faster rate t han surface avadabiht . for 
space radiator*, particularly on reusable vehicles 

A ImpoMtion of contamination control constraint* may pre- 
clude the u-e of water evaporator* or like device* in the vicinity 
of astionomv or earth observation experiments 

In addition, future reduction in spacecraft power eoata and the 
availability of a vapor compression ma< hine in the vehicle, nuch 
as in reentry vehicles which have large cro**-range capability, 
will tend to make refrigeration of -pa< e< raft more attractive 
Refrigeration sy-tems become more attractive in spacecraft a* the 
de-ired coolant outlet temperature approaches the equivalent ra- 
diator -ink temperature, in which case the required space radiator 
area i* dramatically increased Tsuallv the desired coolant tem- 
perature i- 4H deg K, which i" required in the dehum idih cation of 
cabin air, and possibly m the TO, removal process If a regen- 
erable life support system is used, the life support system proceaa 
heat load requiring cooling to 40 deg F 27H deg K i* small eoro- 
pa.ed to the total heat load In addition, much of the electronic 
equipment w hi* h is cooled to 70 Ml deg V 2**4-3HM deg K) could 
l»e designed for higher temperature operation Thu- the vehicle 
head load i ould t** partitioned into low temperature 40 deg F l 
i JT »s deg K an«i high temperature lm> deg F [dll deg K] 
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section^, and the high temperature section would contiun most 
of the load Foi the space shuttle orbiter, the split could he 20 
!M) percent low tempei attire and 70-SO percent high tcnq>erature 
Of rmirM*, a low tempeiuture sink would almost always he 
available at some tunes for a spacecraft in earth orbit, and at 
thoee times the radiator could operate a> a conventional pumped 
fluid system. 

Hybrid Systia Dtveiopntit Prsirai 

Kig. 1 shows a schematic of the hvbud conventional ladtu- 
tor refrigeration cycle heat rejection system, a> configured for 
the space shuttle orbitei. I'nique development problems asso- 
ciated with tfu> device include fluid inventory control and 
operating mode transfer procedures, in addition to the uMial ones 
associated with two-phaae flow in zero-g 

The mode transfer operation is crucial to the success of the 
system During transfer from the liquid single-phase mode to 
the vapor compression mode, the compressor is user! to force 
liquid out of the system into the accumulator This ojieratton 
requires that the eompret*or pump liquid for a short duration, 


thii' the I'ompieosor iua\ need a variable sper.i drive mi that no 
d:uu ig»* will <kiui duimg the modi 1 tian^fer o^iei-tiion The 
lelativcU -in ill pi a linn ot the heat lo;wl expected to be refrig- 
erated dess than . »m.im M i H'T t ' hr }s 7tf k\V| foj the shuttle or 
biter), couple! with the nen*s-it\ for high efheieney m » > re«luce 
power requirements, mdn at*- the use of a jnisitive displacement 
compressor An efficient «entnfug:il niiuprts«or designed for 
the required heat load would have toofieinie at minting speed* 
of above 1U0,00H tpm 1 1 j .• and if i- nnlikeU that it could survive 
the liquid pumping operation without damage In contrast, 
|M)s|tive displacement minpit -sot- have lit lie difficulty with liquid 
slugging, espenallx :ti Inw q>eed < uididate compressors inrlude 
bl i reciprocal iiig, with reed valves, i h , rotary mi vane i, with 
reed valves, in heh-rotot The requirement for operation in 
/.ero-gravit v will also afTei t the < umpnsMn selection tJenerallv, 
compressor* used in terrestrial operation* have gravity actuated 
lubrication s\ stems However, the rotary ami heli-nuor coni- 
pressois (an be designed to obtain lubrication from the flowing 
refrigerant, either fiom oil si t , ; tended in the refrigerant, or from 
the tefi igeiant itself Theiefore, these compressor' do not re- 
quire gravit\ initialed lubrication s\^tenis, and, in fact, such 
compressors have !>een User! in high [a-rformame aircraft which 
sustain advei'C gravity envitonments foi jxuiod- of do 40 sec 

For feasibility toting to l>e conducted in a terrestrial environ- 
meiit compressor selection was made on the basis of low co.st and 
availability, it was thought that some » ompressors might l>e 
damaged bv liquid 'lugging as sy tern operating technique* were 
|mt erted The r ompressoi selected was a Tectim*eh M<»del 
lPH), whn h |s normally Used in automobile vapor cycle systems 
It is a one-cylinder re< iprocatmg compressor with reed valves 
The coinpressitj ojeuated succes,sfuli\ through jieriudic intervals 
during a year of testing, including -ome lOO operating mode 
transfer operations, without degradation in overall system j>er- 
formance 

The refrigerant selection was related to the decision to tise a 
mass produced u . , low net t compressor, and refrigerant 21 
i H-12 ) was selected on the basis of lowet vapor pte-Mire than the 
other refrigerant t K-22 < widely used tn small |H>s.tive di*place- 
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ment compressor* R-12 has not been u-ed in conventional 
spacecraft radiator system*. principally berniM* of a relatively 
low -pec i fic heat ami a lelaiivelv high vapoi pressure m- com- 
pared to RSHOa, R-21, etc j In modern spacecraft coolant 
selection, freezing temperature i> one of the nuM significant 
parameters, anti R-12 hit- an excellent freezing point of — 2.32 de^ 
F (ll.i..) deg K^. Therefore, it could be u>->ed a* the flight ie- 
frigerant coolant R-21 has received wide attention as the 
coolant for future radiator s% tern- J2, d], and since rotary com- 
pressors are available for R-21. it could also be used as the re- 
frigerant in a flight \riMonnf the hybrid conventional radiator 
refrigeration system R-21 n an excellent refrigerant and has a 
12 percent higher theoretical coefficient of performance (COP) 
than R-12 | 4 ). 

The accumulator also repi events a majo- development item for 
the hybnd radiator refrigeration system It must be capable 
of maintaining system pressure at a different level for each 
operating mode, and it must be capable of supplying or storing 
the refrigerant inventory difference between the two operating 
modes. A bladder-type tank pressurized from an external nitro- 
gen supply was selected for the accumulator A rejected NASA 
Smut launch Vehicle attitude control system propellant tank 
and silastic bladder were initially used in the system testing 
The bladder failed due to t. flaw (which had caused it to be re- 
jected on Scout }, and was not replaced m «u Sequent testing in 
which gravity was used for separation of liquid and the pressiirant 
gas. 

A prototype zero-g evaporator (shown in F.K 2) was used in 
the feasibility tests. If con*isf* of a fube-in-tube heat ex- 
changer, and has swirl fin- in it The evaporator picks up heat 
from a separate fluid loop which flows m the outer tube. This 
device i" a “dr>” ev&pnrntoi, \ , it has continuous flow as op- 

posed to the “flooded " evaporator where the liquid refrigerant 
accumulates in a pool around the coil and boils off on demand 
Refrigeration systems in high performance aircraft normally Use 
“dry" evaporators (’out ml* foi this system were not developed 

during this woik, manual mntiol of all system elements was 
generally used, the exception being that the constant superheat 
expansion valve wa- allowed to contiol evaporator outlet tem- 
perature during some te-t pha-e- 

The bulk of the activit) in thi- work was devoted to develop- 
ment of the condensing tadiutor, and the remainder of this prper 
us devoted to that effort 

Rarfiattf/CMfnur OasifR 

\n existing radiatoi design pj] was used a^ a prehnunaty 
model for the R( ‘ busting anahtical technique- and experi- 
mental data were used a.- a ba-i- foi zero-gravits , iwo-phtt.se 
flow design calculations 

After finalization of the RC design, two R(’ panel- with tul»e 
diameters of 1 Vm [.i [s mm) and (MHj.Vin J2 4 nun , wete 
fabucated and verification testing of the entue *\-tem wa- 
‘ onducted under one-g condition- The R(’ panel with 0 »HC»-m 
2 4 mm) tubes was subsequently tested in a C-blo aircraft 
under simulated zer<eg conditions In addition, flow visualisa- 
tion tests of a model of the RC panel, involving high speed 
cineinai.ifraphy, were conducted under one-g and simulated 
zero-g conditions 

The design goal' for the condenser radiatoi i RC . weie 

a Have the same external dimensions a* the individual shuttle 
orhiter radiator panels fi-ft ■ 1 Mm ■ wide l>\ 12-fi i t.fim in 
length t 

6 Operate at a maximum conden-ing tempeta’uie of I deg 
Y t Xis q deg K 

e Reject up to HTT hi 2*07 k\\ m 4 -u,k t.-uqier.*- 
ture of ti deg F 27d 1 » deg K 

n c ion -taut area parallel flow p 4 --age- ..f 4 b«. n \ s lfl 

. - mm 1 1) [2, 




TU« CROSS S£CTlON 

sural fins in ioth tu«$ 

MAT l 

comer throughout 
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r Have a ftowiate of at lea-t 40 lb hi-tuhe < *> 04 g —tulteiin 
the all-liquid operating node sufficient to provide turbulent 

fti »W' ) 

f Incur a maximum R(‘ piessiue dmp of p*i dCi kl*»i in 
the all-liquid mode at, 1 4 p-i ( 2 s k Pa - in the vapor-r\ « le mode 
g Maintain hvdrodvnamtcallv stable flow in the UC in the 
vapor-cycle mode, such that no va|n>r is introduced into the RC 
outlet manifold 

n Provide positive liquid collection at the outlet heade- 
i Operate successfully with an acceleration of O 1-g along any 
ax is of the i ad i a tor 

Satisfaction of these design goal- was believed ne<e*saiy in 
yield a heat rejection -s*tem which could ea.-ih he integrated 
into current or advanced spaceci aft 
Tww-PNwt* Hww SwWIity. There has been considerable investi- 
gation of condensing two-phase flow stability, in one-g, ami some 
also in zeto-g fields Several t\ pe> of two-pha*e flow can exist m 
tubes, examples of which are shown in Fig l Literature mu- 
veved indicated that annular flow i- the most -table two-phase 
flow regime foi a wetting fluid Theteims ‘wetting and “non- 
wetting’ refer to the a< H<»ti of a liquid on a surface For example, 
liquid mercurv read- to ball-up on alum. . urn, so it i- a ‘'non- 
wetting" fluid , R-IJ irtids to spread out *o it in a wetting 
fluid * HasicalR stabihtx, in the context of RC operation, 
means that the flow regime patterns along the length of the tubes 
hie ron-taii* with time **i|ch behavior tends to prevent large 
inventors changes m t h*- RC, whiifi would affect the operation 
of t he -\ - f em since over all inventors t* fixed for each operating 
mode It al-o tend- to prevent laige pressure flue* uationx, and 
hlialls o-i illation of * he liquid leg of the tuln-s is avotde^i If 
oscillation ociiu* M * ould result in vajeu p t uitg imo tfie liquid 
header, **t plugging of the »uf*e fq\ liquid f>a* kifig »fi * « » tfie vapor 
fieadet Hn adls -(jerking, the con ht tori of -t 4 t, tilts n achieved 
It per t m f»at ions ’he -s-tem 4 je dftfiqied ou* i at het ? hati l*eiug 
amplified 1 fie most -ignifi* an* ? sjh*> * f .innular fl«>w iri-tabih- 
tie- are interfa> iuI m-tibihts **t die vaj*of liquid in'rtbo* 
and film oi KsdoMi^tiam iii-tifidi’s 

l he I U-sh tot m f> u at.nulai * ornlet. -mg fl“W t - -flow it m F ig 4 
'1'he liquid film remain- n I »Msel\ ditf,, i h v- m* te.t-mg s eioi i’ s 
m * fie dite< • joj, . .f fl> > w w fill** ’fie \ t (ci s el< < i * s. > fe< i en-e- 4 s tfie 
vaj.**! In t-> tloW l Ot de If I'* * ipidi . ’f.Jtfi die \Ajx-i fl • » w 

ii*i I S * n 1 • i it ' ’ fu \ ■! ■ e t n* i Id iU 1 , u ’ ei t * i* ftei e 1 1 1 

* he s 4 {m.< i- < * i*u. • t V’ *f tl p«.,r • *fie liq u*t vel" i* t i- 

q i.te i 1 1 W if,- « ’ 1 w q>,» ’f* * i . ’ 1 1 e fl* w * ft iliftel 
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critical diameter for R-12 mm m function of temperature. Kor an 
operating temperature of 150 deg F (33M.9 deg K) the critical 
diameter in a one-g field m 0.045-in. (1.14 mm), while in a *ero-g 
field it ia infinite There are numerous types of film or hydrody- 
namic instability, such as : 

1 Tollmien-iSchlichting, which is a result of viscosity amplifi- 
cation of disturbances occurring when the liquid film undergo** 
transition from laminar to turbulent flow (5, 6j 

2 Kelvin- Helmholtz, which results from interaction between 
the vapor and liquid films in the form of aerodynamic premure 
which tends to lift the liquid (7] 

3 Gravity runback, which occurs when an acceleration force 
along the axis of the tube oppoaite to the direction of flow causes 
the liquid film to flow back against the vapor flow 

The biggest factor leading to liquid film bridging and develop- 
ment of slug flow appears to be the Kelvin-Helmholtz instability 
17}. Miles (9| has developed a two-phase flow stability criterion 
which may be expremed m the form of *h* liquid *’m Reynolds 
number and the V' ber number at the liquid surface Basing the 
liquid Reynolds nuinoer on free surface velocity and film depth, 
Miles determined that for large liquid Reynolds numbers, i.e., 
Re* -* * , a sufficient condition for stability is thai the Weber 
number be less thn;. d As a lower bin d, he found that there 
exists a critical liquid Reynolds number below which the flow 
will always be stable. His calculations indicate that the value of 
this Reynolds number w approximately 200 These two stability 
criteria are calculated as follows 



4 C— d—wgn wMi • tMrfe 4aw 


This liquid velocity decline is accompanied by a static pressure 
increase, due to the change in momentum Interfacial or Taylor 
instability can occur at ihte point when gravity causes the inter- 
face to collapse, allowing vapor to proceed beyond the interface 
location This type of instability is obviously a function of the 
tube diameter and the local acceleration field Demngton (MJ 
has defined the maximum stable tub* d imeter a* 

!A < 1 *15 I n <h 

L'p p.*J 

Surface tension <•> ami defi'it) p are properties <4 the par- 
ticular fluid, and are a function of tem(ieraiurr Fig S shows the 


Xw * < 3 for stability (2) 

Rei - PL “* * < 20 0 fur stability 13) 

Mi 

It can be seen from equation* < 2 » and that the criteria for 
hydrodynamic stability during angular flow mat t>e determined 
by “adulating liquid free surface velocity and film thickness. 
Since these parameters change over the entire length of the 
condenser tube, the Weber and liquid him Reynolds number 
mu»t be evaluated over the entire tube length 

This t riterton wv# applied to the present design for a flowrate 
(4 20 lb hr f2 52 g ». i, which represents th • flow at (tie maximum 
heat toad in the vapor-cveic mode, with the results shown in 
Fig 6 Thew* r.*»ultn indicate that a !J s in <9 .Vi min i I 1) tube 
is required for stability over al* liquid-film Reynolds numbers 
expected for a 12-ft d 66 etn i tube The maximum liqui** 
Reynolds numlter approaches 2hfJH t so increasingly smaller 
diameter tut>*s would tend to pi since m-tabilit v prior to com- 
plete condensation These results tend ti -how that stability 
< annot l»e a* hieved with the desired tid** diameter of I Vm fils 
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mm). This small diameter i? desirable because it; fa) is the op- 
timum diameter for conventional radiator operation; (6) favors 
system operating stability because of a smaller fluid inventory 
in the condenser; and (c) requires a smaller accumulator for the 
operating mode changeover, since the fluid inventory is smaller. 
For these reasons, the flow regime for the 1/8-in. (3.18 mm) tube 
was explored further. 

The approach used to predict the flow regimes expected within 
the condenser tubes was based upon the experiments of Baker 
[10] . His results are presented in the form of a flow regime map, 
shown in Fig. 6, and are based upon adiabatic two-component, 
two-phase flow. Since two-component, two-phase flow results 
have frequently proved to be applicable to single component, 
two-phase flow systems, it was considered that use of Baker’s 
flow-regime map was justifiable. This method is based on the 
parameters; 


[U)Gs)r 

(4) 

(?)[-(“)? 

(5) 


and plots of (7,/X versus GlX^/(?,. Fig. 7 shows a Baker Chart 
with various flowrates for a 1/8-in. (3.18 mm) tube plotted on it. 
in the condensing flow case G * = /(C?l), and if it is assumed that 
quality varies linearly along the condensing length, then 
since 

Gr = <?» + Gl (6) 



Fl|. 6 Hydrodynamic stability criteria 
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X = 


G v 

Gt 


(7) 


and 


X & l/Lc (8) 

it is possible to solve for Gl ! G v as a function of L c by combining 
equations (6), (7), and (8) 


Gi _ L e - l 
G v i 


(9) 


This then defines L c as a function of G v /\, This function is 
shown on Fig. 7, and it indicates that for a 1/8-in. (3.18 mm) 
diameter tube, a transition to slug flow will occur at about 1 /2 
of the condensing length, and that plug flow will occur at about 
0.9 of the condensing length. 

Thus, these analytical results indicate that annular flow will 
not exist down the entire length of the condensing section. 
Possibly as a result of simile* analyses, many authors [6, 8, 16, 
and 17] have suggested that tapered tubes are necessary to pro- 
duce stable flow in zero-g condensation. They reason that re- 
ducing tube diameter in the direction of the flow increases the 
vapor velocity, thus helping prevent bridging of liquid across the 
tube. It has also been suggested [6] that the use of orifices in 
each tube at the liquid manifold would produce a more stable 
flow. These approaches were demonstrated to have merit by 
Soliman [6] . However, in the present work the tapered tube was 
judged too difficult to produce, and introduction of orifices was 
found to increase the pressure drop appreciably 

It is possible for stable slug flow to exist, of course [11 ] . This 
seems to be a particularly reasonable expectation for small 
diameter tubes, where capillary action is significant. Work done 
on slug flow by Suo p2, 13] indicates that stable transition to 
slug flow occurs in capillary tubes, and based on criteria given in 
that work, the flow in low quality regions of the tubes presently 
under consideration should be in a region of dominant capillary 
action, particularly when under zero-gravity conditions. On 
this basis it was decided to use straight tubes with no orifices for 
the initial testing, and to resort to the use of orifices only if it 
were demonstrated that stable flow could not be achieved without 
them. 

Profttur* Drop Computation*, \ccurate prediction of two-phase 
flow pressure drop in tubes is extremely difficult, and requires 
knowledge of the flow regime along the entire tube. However, 
the most widely used technique for predicting two-phase flow in 
tubes is that of Marti nelli [14], which assumes that the pressure 
drop of the separate liquid and gas phases are equal, and makes 
no assumption regarding the specific flow configuration. Mar- 
tinelli's predictions assume that the flow is steady and that an 
annular, dispersed flow pattern exists [15]. A computer routine 
which calculates a steady-state flow and thermal balance on a 
single RC tube was developed for use in RC panel design. This 
routine is described briefly in Appendix A. It uses a modified 
version of the Martinelli approach (described in Appendix B) to 
compute the two-phase flow pressure drop in the tube. The 
liquid film thickness is computed from the Martinelli correlation 
for annular flow as described in Appendix B. Two-phase flow 
stability parameters are also computed in this routine. Total 
tube pressure drop is calculated by use of Darcy equation for the 
single phase regions. 

Results of the pressure drop analysis are shown in Fig. 8. 
These results indicate that the 1 S-in. (3. IS mm ) and larger tubes 
yield a pressure drop which is within the required range of less 
than 4 psi (28 kPa ). 

Hoot Tran •for. As described in Appendix A, the computer rou- 
tine calculated steady-^tate heat transfer for '•mail elements of 
the tube in sequence along the entire tube length The familiar 
1 >ittu^-Boeltei equation wa* u**ed for determination of heat trans- 
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fer coefficients in the turbulent single-phase flow regions of the 
tube (Re > 2000) 

h = 0.023 Re° 'P f » « (10) 

and 

h - 3.66^ (11) 

was used for laminar flow (Re < 2000) 

Several techniques are available for calculating the heat trans- 
fer coefficient in condensing flow [18-20]. However good agree- 
ment was not found among these approaches. In this case, 
the condensing heat transfer (hA ) term is one or two orders of 
magnitude greater than the radiation (UA) term, so great ac- 
curacy in the condensing heat .ransfer coefficient is not required. 
The method of Akers and Rosson [9] was selected for use in this 
work, primarily because it was available in coded form from 
another computer routine. 

It was planned to design for a minimum area radiator, rather 
than a minimum weight radiator, so radiator fin efficiencies of 
0.9 [24] were considered. For the 0.063-in. (1.6 mm) thick 
aluminum plate selected for the radiator (on the basis of ease of 
fabrication) this resulted in an 8-in. (20.3 mm) tube spacing. 
For the 72-in. (1.83m) wide panel, this results in 9 parallel flow 
paths. The radiators were coated with 3M Black Velvet Paint, 
which was assumed to have an emittance, <, of 0.9. 

Results from the routine are given in Fig. 9, where required 
condenser tube length is plotted against inlet saturation pressure 
(and temperature) with inlet temperature as a parameter for 
various flowrates. These results show that a 12-ft (3.66 m) tube 
will provide complete condensation for a flowrate of 20 Ib/hr-tube 
(2.52 g/s-tube) in a 1/8-in. (3.18 mm) tube with an inlet pressure 
of 220 psi (1517 kPa) and an inlet temperature of 250 deg F 
(394.4 deg K). The results also show that the required tube 
length increases with increasing inlc* temperature, and decreases 
with increasing inlet pressure (that is, with increasing saturation 
temperature). 

Mietor/CeMwtMr Panel Datlfft. Based on the foregoing analy- 
sis the radiator configuration shown in Fig. 10 was evolved. It 
was decided to test two RC panels with different diameter tubes. 


Both RC panels were made from 0.06^-in. (1.6 mm) thick 5052 A1 
sheet, with extruded 6063 A1 tubes. The tube diameters were 
t/8-in. (3.18 mm) I.D. and 0.095-in. (2.4 mm) I.D. of configure 
tions shown in Fig. 10. A 1-in. (2.54 cm) I.D. gas header and a 
1/2-in. (1.27 cm) I.D. liquid header were used on both radiators. 

Test PreiriR 

Three series of tests have been conducted and a fourth is 
planned. The first tests involved the complete hybrid radiator/ 
refrigeration system in a space environment simulator (SES). 

Hybrid Radiator R d il mdlw fyttore Ttottof. Fig. 11 shows the 
test set-up with the two radiator panels mounted back-to-back in 
the SES. The panels were tested sequentially and each was 
tested in the radiator and the vapor cycle modes. There were 
also numerous mode transfer operations made during the course 
of the test. 

Fig. 12 compares measured and predicted RC panel tempera- 
tures along the length of the panel for two flowrates. These re- 
sults are typical and demonstrate the general good agreement be- 
tween predicted and measured results that was obtained. The 
knee in these curves indicates the point where condensation is 
complete. The point of ocmplete condensation is about the 
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Fig. 11 Rodlotor/conduntor ponolt mounted in space environment 
simulator 


same in both predictions and data; however, the data do not 
show as abrupt a change in temperature as the analysis. The 
analysis is undoubtedly weak in the transition region, since it 
assumes annular flow throughout the entire two-phase flow 
region. Predictions of flow regime, however, indicate that slug 
flow develops, but there is no suitable heat transfer correlation 
available for this region. 

Pressure drop predictions are compared with test data in 
Table 1 These results are especially good since prediction of 
two-phase flow pressure drop is usually only accurate to about 60 
percent [15]. In this case not only was the two-phase flow pres- 
sure drop predicted, but also the two-phase flow length was in- 
cluded in the prediction. The predictions are consistently high, 
and this may be due to the flow regime model used in the analysis 
which considered annulai flow to the transition point. If this 
is true, then the liquid bubble type flow may have a lower pres- 
sure drop than annular flow. In any event, the Martinelli 
assumption of equal pressure drops for each phase is not true for 
plug or large bubble type flow. 

Another element of the test involved elevation and depression 
of the outlet of the panel to simulate 0.1-g acceleration along the 
axis of the radiator. No flow instability was evident from this 
test. Against the adverse gravity, the condenser pressure in- 
creased and the condensing length decreased, corresponding to a 
longer liquid leg, while in the downward flow case, the condenser 
pressure decreased and the liquid leg was shorter. 


DISTANCE FROM 
INLET - M 



Fig. 12 Comparison of prodlctol and moosurod rodlator/«ond»n»ur 
oanol Ism potato r* profito 


Flow Visualisation and Zoro-g Totting. The SES tests yielded 
good system performance data; but it was difficult to obtain in- 
sight iUo the nature of the two-phase condensing flow with the 
resulting temperature and pressure drop data. The most direct 
method of obtaining this insight is by visualization of the flow, 
and use of very high speed photography to slow the action to a 
point where it can be observed. Visualization allows positive 
identification of the flow regime encountered. It also clearly 
chows if there are any gross instabilities, such as oscillation of the 
liquid leg or similar phenomena. In addition, it provides data 
on liquid surface wave frequency and amplification which can be 
used to upgrade instability theories. Also, data under zero-g 
conditions were desirable to provide confidence that the RC 
panel would perform as expected. 

For these reasons, a test program involving simulated zero- 
gravity conditions produced by a C-135 aircraft flying trajec- 
tories such a.s shown in Fig. 13 was initiated. This type of simu- 
lation produces a low-acceleration field for 15 to 20 sec, and has 
been used extensively by Namkoong and others in studies of 
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non wetting flow for use in design of spacecraft power plant RC 
panels [21-23], 

The time interval for establishment of any change in condi- 
tions is short, but this should be sufficient for any hydrodynamic 
changes to the flow From Dennington [S], the time for a 
surface in a tube to obtain a new equilibrium shape when it 
moves into zero-gravity is 


e 


L — - 

J 2<rg c 


( 12 ) 


Thus, according to equation (12), the time required for the sur- 
face in a 1/8-in. (3.18 mm) ID tube to move 6-in. (15.24 cm) is 
0.66 sec. 

The time for establishment of new equilibrium thermal condi- 
tions is less well defined, hut the change in thermal conditions 
was believed to be of such small magnitude as to have only a 
second order effect on the flow pattern. 

The flow visualization experiment was designed as shown 
schematically in Fig. 14. R-12 is stored in a tank which can be 

heated to provide the desired system pressure. Liquid R-12 is 
withdrawn from the tank and put through an electric heater 
where it is raised to the desired temperature level. The R-12 
then flows to a gas receiver which damps out any pressure 
fluctuations caused by flashing in the heater. It then flows into 
a heater with three parallel tubes made of fused quartz, down the 
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Fig. 13 Aircraft zero gravity maneuver 



Fig. 14 Test apparatus flaw schematic 



Fig. 15 Flaw passage arrangement 


8-ft (2.44 m) length of the tubes, and then overboard. The cool- 
ing medium is water which is stored in two bladder tanks pres- 
surized by aircraft nitrogen supply. The water flows in parallel 
with the R-12 in a channel comprised of 1/2-in. (1.27 cm) thick 
plexiglass, as shown in Fig. 15. 

The complete experiment package is shown in Fig. 16. The 
high speed cameras are mounted on optical benches which permit 
them to be moved along the axis of the tubes. Two cameras 
take motion pictures of the outside R-12 tubes, while a third 
camera photographs all three tubes from above. The flow 
passages and the optical benches are shock mounted to damp out 
all vibrations above a frequency of 12 Hz. 

As indicated in Fig. 14, the RC panel with 0.095-in. (2.4 mm) 
ID flow passages was also included in the aircraft testing. It 
was instrumented with 27 thermocouples, 4 absolute pressure 
transducers, and 5 pressure drop transducers. Heat was re- 
jected from it to the cabin atmosphere and surroundings. 

The visualization experiment indicates that the two-phase 
flow regimes down the tube are as shown in Fig. 17. There is a 
considerable amount of wave action on the liquid film from the 
earliest observable point. The flow gradually develops large 
waves which eventually bridge the tube and isolate large vapor 
bubbles. These bubbles flow rather steadily down the tube, 
while growing steadily smaller until they are finally extinguished. 
This same flow pattern prevails over most of the range tasted, 
that is, from 5 Ib/hr-tube (0.63 g/s-tube) U) 20 lb/hr-tube (2.52 
g/s-tube). At 5 lb/hr-tube (0.63 g/s-tube) it appears that the 
flow is initially stratified, rather than being annular. 

The flow is steady in appearance, and there is no sharp 
boundary between the condensing region and the pure liquid 
region. Only minor oscillations of the liquid leg were obeerved. 
The fluid inventory in the tube gives every indication of being 
remarkably constant. 

The pressure instrumentation in the tests conducted so far has 
been quite erratic, but it appears that there are no significant 
pressure fluctuations. During the portion of this test where the 
outlet end of the tubes were raised to represent an adverse pres- 
sure gradient, the liquid leg length increased by 100 percent. 
However, the pressure instrumentation was not operative at this 
time, so pressure changes could not be recorded. High-speed 
photography of the flow indicates it is not substantially affected 
by the adverse gravity gradient. 
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Fig. 18 shows a sequence of photographs of the flow at different 
stations along the tube; the photographs were not taken at the 
same time, but are at essentially the same conditions. This 
sequence illustrates typical flow conditions observed in one-g, 
and in zero-g. 

In zero-g the flow appears to be more annular in nature than in 
one-g, with a noticeably smaller liquid film on the bottom of the 
tube. Bubbles in the flow tend to be centered in the tube, 
rather than riding along the top of the tube as they do in 1-g. 
Also, the length required for condensation appears somewhat 
shorter than in the one-g condition. Additionally the flow has 
the same stable character in zero-g as was observed in one-g. 

Condusious 

The primary conclusions reached in this study were that a 
hydrodynamically stable RC can be designed for operation with a 
spacecraft refrigeration system, and that a hybrid conventional 
radiator /vapor cycle heat rejection system is feasible. 
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APPENDIX A 

Csntfmini Radiator Analysis Computer Routine 

A computer routine was written for the purpose of analyzing 
condensing two-phase flow in a space radiator. The routine 
provides a steady-slate solution, which is generated by doing 
successive steady-slate heat transfer analyses on small clement* 


along the length of a single tube of the radiator starting »t the 
inlet as shown in Fig. 19. It is assumed that there is an adiabatic 
mid-point in the fin between the tubes, and that all tubes are in a 
similar environment. This appendix describes the computer 
routine very briefly. 

Input Data 

O Inlet fluid flow conditions 
— Temperature 

— Fieaeuie 

— Flowrate 

O Environmental conditions 
— Sink temperature, T , 

— Fin surface emittance, c 

O Radiator geometry 

— Internal flow passage diameter, D 
— Fin thickness, t 
— Fin width, W 

— Tube length, L (The routine will continue to incre- 
ment length until condensation is complete if a 
zero is input for length) 

O Computation parameter 

— Element length, &L (Normally 0.01 ft [3 mm) was 
used) 

— Output frequency along tube length (usually 6-in. 
[15.24 cm] or 1 ft [30.48 cmj was used) 

lotaruil Computations 

The routine was specifically set up for refrigerant 12, and it 
calculated all required fluid properties such as saturation pres- 
sure, viscosity, surface tension, etc., at the required temperature. 

The sequence of computations was as follows: 

(а) Physical properties of fluid based on element inlet con- 
ditions. 

(б) Heat transfer coefficient baaed on the method of Akers 
and Rosson [19]. It was considered to be reasonable to start 
with the heat transfer coefficient and to not iterate by recom- 
puting heat transfer coefficient after liquid film thickness is 
established, because the fin radiation heat transfer dominates the 
heat transfer process. Large changes in the fluid heat transfer 
coefficient have little effect on the overall heat transfer rate. 

(c) Heat transfer from the fin, and steady state temperatures 
are computed. 





(d) If the fluid temperature drop puts it below the saturation 
temperature, the fluid temperature is set at the saturation tem- 
perature, and mass transfer from the vapor to the liquid state is 
calculated. 

(c) Pressure drop is computed using the Darcy equation in 
single-phase flow, or using the modified Martinelli approach de- 
scribed in Appendix B for two-phase flow. Liquid film thickness 
in two-phase flow is calculated as described in Appendix B. 

(/) Liquid film Reynolds number and Weber number are 
computed. 

(p) Steps a through f are repeated for each element of the 
tube in sequence. 

(h) The total radiator pret. re drop is computed. 

Output Data 

Results are output at specified locations along the tube length 
as follows : 

Fluid temperature 
Tube temperature 
Liquid film thickness 
Gas Reynolds number 
Liquid Reynolds number 
Liquid Weber number 
Gas flowrate 
Liquid flowrate 
Incremental pressure drop 


The correlation assumes that an annular dispersed flow pattern 
exist", with no knowledge of other flow configuration details 
being required 

The key to the approach is the identification of a flow modulus, 
a, which accounts for the geometry of the flow, as well as relative 
motion or slip effects between the liquid film and the gas core. 
The flow modulus can be defined as 


a 



(13) 


However, the quantities A L and Dl cannot be directly deter- 
mined, in general. 

In this study, a stability evaluation of the two-phase flow was 
desired, and this required a knowledge of the liquid film thick- 
ness. The most stable flow patw.n, and thus the most desirable, 
was annular flow. If the annular flow were stable, then, in a 
zero-gravity environment, it is possible to express a in terms of 
the liquid film thickness, 5. This is because 5 is the same all 
around the tube periphery, which it would not be in one-g, or if 
the flow were less stable. 

The method for calculating pressure drop and film thickness is 
developed as follows from reference [14] : 


/AP\ /A P\ / •• 

\A L)t-, " \A L), \dJ 

(U) 

For turbulent gas and turbulent liquid 
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(15) 

For viscous liquid and turb ulent gas 
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(16) 

where 

‘-teresr® 

07) 

Cf, Dm If L 

A., - -*<!*,>—« = W 

C, pl M* "» 

08) 


In addition, the following are output for the tubes: 

Outlet temperature 
Outlet pressure 
Total pressure drop 
Condensation length 

APPENDIX B 

Molilil Martiiilli Mithai far Calcalatiaa al Twa Pkau 

Flow Prassura Drop and liquid- Fils Thickaass 

Prediction of pressure drop hr two-phase condensing flow in 
a space radiator requires a suitable analytical model of the flow. 
When two phases flow concurrently in a tube they can assume 
man) complex flow patterns, as shown on Fig. 3. To accurately 
predict two-phase flow pressure drop, mathematical models 
ba^ed on a particular flow paucrn are required, since no single 
correlation or mathematical model can p ove/ all possible flow 
pattern^. Due to the complexities involved in » purely mathe- 
matical analysis, experimentation is used to obtain solutions. 
A well known correlation was developed by Martinelli [141, in 
which the prevmre drop of each phase was considered equal, and 
correlating parameter** were derived from relations governing the 
pressure diop of each of the pha^e* a* if it were flowing alone. 


In this study the gas was always turbulent (i.e. f Re > 2000), and 
the liquid transition Reynolds number was 1000. 

For the assumed perfect annular flow 


d l - 24 

(19) 

D, - D, - 24 

(20) 

D, - 4 
4 

(21) 

Dr c* + 1 

— — 

(22) 

Di a - 1 


Substituting equation (22) into equation (15) 

(<pr~jY ’ - <* + a* ‘ •'«' •)■ • ( 23 ) 


which reduces to 


Am 


6.35 

(a - 1 )' mi 


(24) 
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By equating equations (17) and (24), it is possible to solve for a 
directly. Then equation (21) can be solved for 6, and equation 
(15) can be solved for ju*, which permits computation of two 
phase pressure drop with equation (14). 

A similar procedure can be used on the viscous-turbulent flow 
to yield 
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Abstract 


High-speed photographs wre taken of the flow condensation of refrigerant-12 in small diameter quartz 
tubes (3 mm). Photographs obtained under 1-g conditions permitted comparison of the flow regimes with 
the Baker chart predictions. Qualitative comparisons with other photographs taken under zero-g conditions 
appear to indicate no significant reduction of heat transfer under zoro-g conditions. 


INTRODUCTION 

The present study 1 was undertaken as part of a 
broader program to develop a hybrid spacecraft heat 
rejection system, i.e. one in which the radiator 
could function either (1) conventionally with 
single-phase liquid flow through the radiator tubes 
or (2) as a condenser in a vapor compression refrig- 
eration cycle, capable of accomodating two-phase 
condensing flow. The impetus for the development 
^f such a nytriu system originates principally 
from tn« trend toward internal power requirements 
of spacr vehicles growing at a faster rate than the 
surface availability of space radiators. Feasi- 
bility studies 2 indicated that a vapor-compression 
refrigeration system could reduce radiator surface 
requirements significantly in some apt’ ' ications . 

While unique developmental and system operational 
problems are associated with the successful opera- 
tion of such a hybrid system [l], additional 
possible problems related to two-phase condensing 
flow processes (within the radiator /condenser 
(r/c) panel) in a zero-gravity environment must 
alec be accurately evaluated and predicted. The 
broad objectives of the present study were (1) to 
experimentally det -mine operating characteristics 
of a full-scale mux ti -channel radiator/condenser 
panel, with refrigerant-12 (R -12) as the working 
fluid, in 1- and Q-g environments, (2) to establish 
and/or modify parameters and correlations capsble 
of predicting heat transfer, pressure drop, and 
stability chat -cteristics of such systems in 0-g 
environments. 

In order to help achieve these objectives photo- 
graphs of the two-phssc flow condensation phenomena 
occurring in the condenser tubes under both 1- and 
0-g conditions were desired. This information was 
obtained by high speed photographs of R-12 conden- 
sing in s system of three quartz tubes connected In 
parallel end modelled thermally and hydrodynamically 
after the r/c panel. 

The present paper concerns itself primarily with 
the photographic measurement end analysis of con- 
densation phenomena in e 1-g field. Alec, some 
qualitative comparisons are made with photographs 
taken in 0-g. Quantitative comparisons with 0-g 
behavior are to follow as more 0-g photographic 
data become available. 


If it is essential to predict flow regimes, so that 
heat transfer, pressure drop and stability charac- 
teristics may be more accurately predicted, it must 
first be determined whether existing flow regime 
asps, such as the Baker chart [2], accurately por- 
tray the flow regimes occurring in the given system 
over the range of operating conditions. Flow 
regime maps are available in the literature for 
both horizontal and vertical flows [2-3]. Gener- 
ally these maps have been developed for boiling 
systems or for two-phase, two-component systems 
(e.g. sir-water), which have often been obtained ip 
rather large diameter systems. Whether or not the 
flow regimes predicted by such maps have any 
validity under 0-g conditions should be eat 'GHshed 
by the present tests. Appropriate modifications to 
such maps, together with the related predictions of 
system performance, may be required. 


PREVIOUS ZERO-G INVESTIGATIONS 

One of the earliest investigst ions of condensing 
processes in zero-gravity conditions was conducted 
by Feldmanls [4], who made temperature, pressure, 
and photographic measurements of forced convection 
condensation of distilled water in s single, 
tapered, tube, with a 1.74 cm inlet i.d. He re- 
ported an increase in condensate temperature during 
low gravity conditions that may have been caused by 
reducsd condensation heat transfer coefficients, 
possibly attributable to the absence of gravity 
forces acting to remove the condensate film from 
the cooled surface within the tube. 

Two studies of significance, investigating the flow 
condensation of s non-wetting fluid (mercury) in 
tubes, were conducted by Albers and Macosko [3] end 
Namkoong, et ml. [6]; both studies employed s Navy 
AJ-2 bomber flying a Keplerlan trajectory to attein 
the 0-g environment. Tube sices ranged from 0.686 
to 1.245 cm i.d. Namkoong, et el. [6] reported a 
concentration of drops on the tube bottom and a 
shallow sloping interface under 1-g conditions, 
while In a 0-g envlroomant the drops were uniformly 
distributed end the Interface wee vertical. Albers 
and Macoako reported that (1) measured overall 
static pressure drops at two mass flow rates were 
essentially the seme for 1- and 0-g conditions, and 
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(2) the Lockhart-Martinelli parameter ^ and the 

two-phase flow pressure gradient were affected 
negligibly by gravity, as also were the fog-flow 
correlations. 

In the study related sost closely to the present 
one, Sollnan and Berenson [7] studied flow conden- 
sation of Treon-113 in a transparent multi-tube 
condenser. Both tapered and constant-diameter 
tubes were studied, the latter containing twisted 
tape and wire in some cases. Tube diameters ranged 
from 0.335 to 0.595 cm i.d. The gravitational ef- 
fects studied, however, wer^ those attained from 
operating the condenser in positions ranging from 
a 4-90-degree position (condensing downwards) to a 
-90-degree position (condensing upwards). 


EXPERIMENTAL EQUIPMENT 

Three test articles were iabrlcated — two full-scale 
radiator /condenser (r/c) panels, (Figure 1) and a 
system of three parallel-connected fused quartz 
tubes (Figure 2) modelled thermally and hydro- 
dynamlcally after r/c panel mmber 2. The r/c 
panel number 2 was intended co provide quantitative 
engineering data in 1- and 0-g conditions, while 
the transparent quartz-tube fixture was Intended to 
permit high-speed photographic analysis of the two- 
phase condensation phenomena, also under 1- and 0-g 
conditions. Testing was performed initially on the 
ground, with both horizontal and inclined flow 
passage orientations, and additionally in an air- 
craft flying 0-g trajectories. 

The transparent quartz tube system is shown in the 
system schematic, Figure 3. The condensing fluid, 
R-12, is stored in a 75.7 liter tank that can be 
heated to provide the desired system operating 
pressure (0.0984 to 0*246 kgf/om 2 ) by an lmmerslon- 
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Fig. 1 RADIATOR /CONDENSER PANEL DESIGN 



Fig. 2 FLOW VISUALIZATION APPARATUS 



Fig. 3 TEST APPARATUS FLOW SCHEMATIC 

type heater. The liquid R-12 withdrawn from the 
tank bottom passes through an electrical heating 
system where the R-12 is superheated to the desired 
degree. The heater is designed with a tortuous 
pathway systos so that centrifugal force* maintain 
the liquid refrigerant in contact with the heater 
walls during weightless intervals. The R-12 then 
flows to a gas receiver so that any pressure fluc- 
tuations caused by flashing in the boiler may be 
damped out. 

The superheated R-12 vapor flow is then directed to 
either the r/c panel or to the transparent tube 
fixture. In either case the exiting liquid R-12 
passes through a flow metering valve which main- 
tains the flowrate at the desired level 2. 3-9.1 
kgm/hr. The condensed R-12 is then either dumped 
overboard or to the atmosphere, depending upon 
whether flight tests (0-g) or ground tests sre 
being conducted. 

The quartz tube fixture ia cooled by water (Figures 
2, 3) stored in two 53-liter tanks in a closed loop 
parallel flow system (though in some initial teats 
an open coolant system was employed). The water 
temperature in the tanka ranged from 15.6-26. 7°C, 
while the nominal flowrate vat 45.4 kgm(m)/hr. 

The maximum temperature rise of the water was 17*C. 
The inlet R-12 temperature ranged from 47.2 to 
86.7°C while pressures ranged from 0.110 to 0.179 

Instrumentation on the transparent tube fixture 
consisted of inlet and outlet temperatures (immer- 
sion thermocouples) end pressures of the R-12, 
together with AP measurements across the quarts 
tube assabiy. The inlet and exit temperature of 
the water system were alto monitored by issicrslon 
thermocouples. Four surface temperature measure- 
ments of the plexlglas flow channel wei* also made. 
Flowrate measurements of both systems were made 
with turbine flow meters, while e rotameter was 
employed In the water system for 1-g teste. 

The photographic observation were made with two 
high speed caaurss (operated at 2000-7000 
frames /second) and a low speed camera (400 
frames /second ) . The low speed camera was mounted 
so as to view all three tubes from above, per- 
mitting observation phenomena within tha parallel 
system. All three cameras were mounted on optical 
benches that permitted rapid repositioning along 
the length of the tubes. The high speed cameras 
viewed the transparent fixture from the side, as 
shown in Figures 2 and 4. 
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Fig. 4 EXPERIMENTAL TESTING ASSEMBLY 


F<g. 6 (c) WAVY -ANNULAR FLOW 


RESULTS 


The types of flow reglnes observed are shown sche- 
matically in Figure 5. This figure has been 
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Fig. 5 SCHEMATIC OF CHARACTERISTIC FLOW REGIMES 


Fig. 6 (d) SLUG FLOW 


developed after examination of numerous high-speed 
photographs taken at periodic intervals along the 
condenser tube's length. The tube inside diameter 
upon which these results are based was 2.62 mm. 

Flow rates of the condensing R-12 ranged from 
2.3 to 9.1 kgm/hr, corresponding to G values of 
42.2 to 168.9 x 10* kgm/m^-hr. Inlet temperatures 
varied from several degrees superheat in some cases 
to saturation temperatures. Over the range of 
variables investigated no substantial differences 
were observed in the qualitative illustration of 
flow patterns depicted in Figure 6. 






Fig. 6 (e) ELONGATED BUBBLY FLOW 
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Fig. 6 CHARACTERISTIC FLOW PATTERNS- 
(a) ANNULAR FLOW 


Fig. 6 (f) ISOLATED BUBBLE FLOW 
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Fig. 6 (b) ANNULAR-STRATIFIED flow 


Figure 6 consists of photographs of the conden- 
sation processes occurring under the conditions of 
a nominal system pressure of 230 psla. Shown in 
sequence are (a) ripply annular flow, with the 
liquid film distributed essentially uniformly along 
the tube wall, (b) annular -strat If led flow, with a 
preponderance of the liquid flow along the tube 
bottom due to gravltat tonal action, (c) wavy 
annular flow. In which the shearing action of the 
higher speed vapor flow Initiates the growth of 
large astplltudc waves, occasionally raachlng the 
top of tha tube, (d) slug flow, in which the regu- 
lar and periodic bridging of the tube diameter 
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results in separated slugs of vapor moving above 
the liquid, (e) elongated bubbly flow, where the 
isolated vapor pockets have become more regularly 
shaped and widely spaced, and (f) isolated bubble 
flow, in which rather small vapor bubbles move 
along the top of the tub’ to complete extinguish- 
ment. 

Evidence of the lack of distinct boundaries between 
various flow regimes is shown in 1 igure 7. Each of 



Fig. 7 VARIETY OF FLOW PATTERNS AT IDENTICAL TEST 
CONDITIONS- (a) ANNULAR -ST RAT I FI KD FLOW 






Fig. 7 (b) WAVY -ANNULAR FLOW 
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Fig. 8 BAKER CHART [2] 


long elongated bubbles to small isolated bubbles). 

Other qualitative resjlts observed primarily in 
connection with differences between 1-g and 0-g 
flow patterns are as follows. As one would expect 
in the absence of gravitational forces, the iondrn* 
sate film in the annular region was uniformly 
distributed about the tube wall. It would be 
expected that the annular regime would persist over 
a longer distance in a 0-g environment. Verifi- 
cation of this expectation and mapping of the 
complete flow regimes in a 0-g environment is to be 
conducted as data from future 0-g flights become 
available. 

Also, as one would expect, It was observed that in 
the bubble flow region the bubbles moved very 
nearly along the tube center, showing no affinity 
for the tube wall. From an overall standpoint, it 
appeared that there was little difference in heat 
transfer characteristics due to 0-g operation since 
there was no substantial changes in the point of 
complete vapor extinguishment for the same flow 
parameters . 

A final qualitative difference between 1-g and 0-g 
flow conditions was the notably less irregular, 
less wavy character of the flow phenomena under 

0- g conditions. In the anrular regime the liquid 
film appeared significantly less wavy. Also, in a 

1- g environment in the annular region, high-veloci- 
tv disturbances appeared to disrupt the liquid film 
appreciably at regular intervals; under 0-g condi- 
tions the periodic high-velocity disturbances also 
persisted, but appeared to disrupt or distort the 
annular film to a much smaller degree. 


I 


L 


the photographs in Figure 7 was taken from the same 
roll of film, which represents less than one second 
of viewing the condensation process. The regimes 
observed at this particular location are (a) annu- 
lar-strat if led flow, (b) same as (a), but with 
complete, but Intermittent, bridging of the tube 
diameter by the waves, (c) slug flow, in which 
distinct separation of vapor slugs by a continu- 
ous liquid column is evident. 

The flow regimes observed a-e compared with those 
predicted bv the Baker chart [j] in Figure 8. The 
coordinate* of 1 igure 8 nre functions of liquid and 
vapor aass flowrate and fluid properties [2,8] 
(viscosity, surface tension, density). The flow 
regimes were categorized Into tour groups: annular, 
annular-st r at it led , slug, and bubbly (ranges from 


DISCUSSION OF RF.SLLTS 

The flow regimes observed, Illustrated in Figure* 

5 and 6, differ substantially from the description* 
ut condensing horizontal flows in tube* presented 
In [4-6,8], Also, as generally described in the 
literature wetting fluids generally form a thin 
annular layer at the heat transfer surface, from 
which layer drops are sometimes broken off and 
entralrsd in the vapor core, resulting In mist- 
annular flow. Very little droplet entrainment was 
evident in the present study, perhaps because of 
relatively low vapor velocities. 

Analytical models of condensation in tubes often 
consider the annular flow regime to persist 
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preponderance of bubbly-flow data point* are a#«n 
to fall along the slug-bubble boundary or slightly 
within the bubbly regim*. In general, the Baker 
chart predictions appear reasonably good. 


throughout the tube length, based upon the suppo- 
sition that annular t low is indeed the dominant 
flow regine. The present results indicate annular 
flow to exist over onlv roughly half of th^ con- 
densing region, after which the slug flow and 
bubble flow regimes exist. 
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In a study of wetting and non-wetting mercury con- 
densat ion [8] f luoroscopic observations indicated 
a fog-flow regime terminated by a rather elongated 
sloping interface at the "point" of complete con- 
densation. Their observation of the non-vetting 
flow Indicated a predominantly fog-flow regime v*th 
some droplets at the walls, terminated by a more or 
less vertical, distinct interface (see also [3,6]). 
These observations led to their employment of a 
fog-fUjv model for predicting pressure drop. As 
seen from r’igure b, not only does a fog-flow regime 
not exist, but neither does a distinct liquid-vapor 
interface at t^t* point of complete condensation 
exist; but rather, elongated vapor bubbles gradu- 
ally become smaller until they eventually disappear. 

Concerning settling of liquid along the tube bottom 
as in Figure h (b) , (c), this regime might be con- 
sidered a combination of stratified and annular 
flow. Regarding heat transfer in this region, 
certainly the accumulation of liquid along the tube 
bottom diminishes the rate of heat transfer, but 
at the same time, continuous draining of rhe liquid 
from the top half of the tube or^'otes heat trans- 
fer in that regiun. Contrasting this picture with 
that existing in Q-g, the essentially uniform film 
thickness existing in 0-g may have a uniformly good 
heat transfer rate along the tube periphery that 
may be equivalent to the "average" heat transfer 
conditions existing in 1-g, described in the pre- 
vious sentence. Perhaps this accounts for the 
observation that no substantial difference in con- 
densation length is discernible for 1-g and 0-g 
conditions for the same heat transfer and flow 
parameters. It should be noted that possible 
differences in the other flow regime also affecting 
the heat transfer characteristics have not been 
investigated due to unavailability of data. This 
will be the subject of future investigation. 

The comparison of measured and predicted flow 
regimes shown in Figure 8 warrants some elaboration. 
The Baker chart is well in accord with physical 
reality in the annular flow region, where at the 
higher system flow rates motion picture photographs 
did give some evidence of the existence of dis- 
persed flow. However, the st rat if led-annular flow 
regime observed has no counterpart *n "he Baker 
categorization. It the st rat if led-annular regime 
may be classified simply as a special case of 
annular flow, then the Baker chart is in reasonable 
accord with reality. 

The onlv three observations of a slug-flow regime 
fall Just outside of the slug-flow boundaries. 

This is not unusual in that the flow regime 
boundaries must be realistically conaidered to be 
somewhat "blurred", rather than discrete. The 


CONCLUSIONS 

The principal conclusions of the present study aza: 

(1) the Bakec chart predicts the condensing flow 
regimes reasonably %;ell; and (2) condensation heat 
transfer is not appreciably affected by a 0-g 
environment for the system described. 


FOOTNOTES 

Conducted under NASA contract NAS9-13410. 

Conducted by Vought Missiles and Space Division, 
LTV, Dallas, Texas, under the auspices of the 
Johnson Space Center, Houston, Texas. 
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